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We investigate the parton-distribution functions (PDFs) for the positively charged pion and kaon 
at a low-renormalization scale ~ 1 GeV. To this end, we employ the gauge- invariant effective chiral 
action from the nonlocal chiral-quark model, resulting in that the vector currents are conserved. 
All the model parameters are determined phenomenologically with the normalization condition for 
PDF and the empirical values for the pseudoscalar-meson weak-decay constants. We consider the 
momentum dependence of the effective quark mass properly within the model calculations. It turns 
out that the leading local contribution provides about 70% of the total strength for PDF, whereas 
the nonlocal one that is newly taken into account in this work for the gauge invariance does the 
rest. High-Q 2 evolution to 27 GeV 2 is performed for valance-quark distribution function (VQDF), 
using the DGALP equation. The moments for the pion and kaon VQDFs are also computed. 
The numerical results are compared with the empirical data and theoretical estimations, and show 
qualitatively agreement with them. 
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I. INTRODUCTION 

It has been well-known that the high-energy scattering processes, such as the inclusive or exclusive production ones, 
are very useful tools to study Quantum ChromoDynamics (QCD) [THTT]. In such processes, the Drell-Yan (DY) one for 
instance, the scattering amplitudes consist of short and long range QCD interactions simultaneously. By the virtue of 
the factorization theorem, the first can be studied via perturbative QCD (pQCD), whereas one is able to investigate 
the second, that signals the nontrivial structures of the hadrons involved, by various nonperturbative approaches. 
Note that those nonperturbative quantities can be defined by the parton-distribution amplitude (PDA), parton- 
distribution function (PDF), fragmentation function (FF), generalized parton-distribution function (GPD), and so 
on, depending on different scattering processes. Among them, PDFs for the pseudoscalar (PS) mesons, (j) = (tt, K), are 
of importance to understand the nonperturbative structure of the mesons, which play a crucial role in the low-energy 
QCD. In Ref. [12 , employing the Nambu-Jona-Lasinio (NJL) model, PDF and valance-quark distribution function 
(VQDF) for the pion and kaon were derived from the forward scattering amplitude of a virtual photon from a pion 
target in the Bjorken limit. The obtained results at a low renormalization scale ~ 1 GeV were evolved to a high Q 2 
by using the Altarelli-Parisi equation. Similarly, in Ref. [13 , the SU(3) NJL model was applied for PDF for the pion, 
kaon and eta mesons with scalar and pseudoscalar couplings. They obtained, in the chiral limit, the simple result 
fcf)(x) = f(/)(x) = 6(x)6(l — x) for the structure functions, satisfying the gauge invariance for the vector current. In 
Ref. [14], the chiral quark model (ChQM) was used for the pion PDF, resulting in that the Pauli-Villars regularization 
scheme is most suitable for both the anomaly structure of QCD and the leading scaling behavior of PDF in the Bjorken 
limit. The pion PDF was computed using the instanton-liquid model, suggesting an analytic expression for a general 
vertex function and satisfying, in a gauge invariant approach [15 . The NJL-jet model was also employed to compute 
PDF as well as FF by the cut diagrams in Ref. [16]. A parameter- free prediction for the ratio uk(x)/u 7V (x) was given 
in Ref. [17] using the rainbow-ladder truncation for the Dyson- Schwinger (DS) equations, reproducing the DY data. 
Extraction of the pion PDF was performed in a next-to- leading order analysis from Fermilab E-615 pionic DY data, 
observing that the high-x dependence is different from that of the leading order analysis, whereas it does not match 
with the perturbative QCD (pQCD) and DS calculations [18]. The lowest three non-trivial moments, calculated by 
the lattice-QCD (LQCD) simulation, were in good agreement with existing data at the physical pion mass [19]. In 
Ref. [20] , the two-flavor Wilson- fermion was used for LQCD simulation, observing smallness of relevant moments in 
comparison with other theoretical models. Using QCD sum rules, the pion PDF was investigated taking into account 
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the nonlocal condensates, paying attention on the bilocal power corrections [21]. In the previous works [22 , 23 , we also 
computed unpolarized FF first then converted it into PDF by using the Drell-Levi-Yan (DLY) relation [24] , employing 
the simplified nonlocal chiral-quark model (NLChQM), although this duality or analytic continuation according to the 
DLY relation is not generally satisfied especially for T-odd PDFs. The numerical results turned out to be compatible 
with available empirical data after the DGALP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) evolution for high-Q 2 
values. 

We note that, although the previous works showed qualitatively reasonable results for PDF and FF [22j [23] . 
from a theoretical point of view, there were some issues to be improved: 1) The vector-current conservation (or 
gauge invariance) of the matrix element for PDF was not taken into account, i.e. only the simple local contribution 
was computed. 2) Moreover, the numerical calculations were simplified to a certain extent by assuming that the 
constituent-quark masses inside the relevant matrix element for PDF are constant, not momentum dependent, except 
for the quark-PS-meson couplings, although this simplification helps to understand the analytic structure of PDF in 
the model and reduces the difficulties much in the numerical calculations. 3) Above all, PDF was converted from FF 
via the DLY, not a direct calculation. Hence, in the present work, we would like to improve all the issues mentioned 
above in the same model. Our strategy in the present work is as follows: 

1. PDF is directly computed by the gauge- invariant effective chiral action (EChA) of NLChQM, satisfying the 
vector-current conservation. This is a similar approach as in Refs. [T5] [25]. 

2. All the momentum dependences in the constituent-quark masses Mj, in which the subscript / stands for the 
quark flavor, are strictly taken into account: Mf = Mf(k 2 ). One can refer for more details in Refs. [26] [27], 

3. Relevant model parameters are determined by satisfying the theoretical condition and experimental information, 
i.e. the PDF normalization condition, and empirical values for the pion and kaon weak-decay constants. 

As already observed in Refs. [25 , it turns out from the numerical results in the present work that the nonlocal 
contribution, which are newly considered here and necessary to preserve the gauge invariance, provides about 30% 
of the total strength for PDF. Similar tendency was also observed in Refs. p~5j [28] in the single-instanton model. 
Moreover, those contributions play the role of broadening PDF. Interestingly, careful treatment of the momentum 
dependences in the effective quark masses does make noticeable changes in the numerical results in comparison to 
those with the simplification, i.e. Mf(k) — >• constant, if we compare them with the previous calculations. Considering 
the normalization condition for PDF and F n = 93.2 MeV as an input, we determine one of the model parameters, 
the constituent-quark mass at zero virtuality, to be Mq ~ 300 MeV, once we choose the model renormalization scale 
as /i « 1 GeV phenomenologically as in Ref. [25 . From this parameter set, the kaon weak decay constant is obtained 
as Fk ~ 121.7 MeV, which is only about 7% larger than its empirical value, 113.4 MeV. From the computed PDF 
for the pion (7r + ) and kaon we parameterize them into a simple analytic form, which was suggested in various 

works for extracting empirical PDF from the experimental data [29] [3D]. Employing the DGLAP equation for high 
Q 2 evolutions for them, we compare our results with empirical data, showing considerable agreement with them. We 
also present the numerical results for the ratio of uk / and the moments of the pion and kaon VQDFs in comparison 
with the data [30 -32 to verify the relevance of the present theoretical work. 

The present work is organized as follows: In Section II, we briefly introduce the present theoretical framework, 
defining PDF and VQDF for the PS mesons. The numerical results and related discussions will be given in Section 
III. The last Section is devoted to summary, conclusion, and future perspectives. 

II. THEORETICAL FRAMEWORK 

In this Section, we explain briefly the definition for PDF for the PS mesons, i.e. pion and kaon. For definiteness, 
we choose only positively charged pion and kaon, and assign them as 7r + ~ ud = tt and K + ~ us = K hereafter for 
simplicity. For other isospin states, it is straightforward to compute the multiplicable factor considering the isospin 
symmetry. In what follows, we want to explain how to define and calculate PDF. In general, information for the 
parton distribution function can be extracted from the (ttotK)-N scattering using the DY process or deeply virtual 
Compton scattering (DVCS) of the PS mesons. The information of the PS-meson PDF can be represented by the 
hadronic tensor W^ v as depicted in the left of Figure [TJ which is so-called the handbag diagram in the forward 
Compton scattering, and it is defined as a function of the Bjorken x by [12j [33] 
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FIG. 1: The left diagram for the hadronic tensor in Eq. for the PS meson with its on-shell momentum p, satisfying 
p 2 — mj. The solid, dash, and wavy lines indicate the pseudoscalar (PS) mesons, quarks, and virtual photon, respectively, x 
stands for the momentum fraction. The thin slashes denote the discontinuity for the imaginary parts. The right diagrams show 
the factorized ones as the hard and soft contributions. The soft part represents PDF with the loop momentum k. The black 
blobs stands for the nonlocal quark-PS-meson vertex, whereas the gray one in the soft diagram for the local operator for PDF 
in Eq. 



where the four momenta p and q stand for those of the PS-meson and virtual photon as shown in the left of Figure [T] 
and v denote the PS-meson mass and (p • g)/m^, respectively. The structure functions in Eq. read 



F 1 (x) = xJ2 e2 f[U( x ) + f^ x )L ^{x) 
/ 



2x 



(2) 



Here, /</>(#) and f^x) indicate PDF for the quark and antiquark, while ef denotes the electric charge for the quark. 
Thus, by measuring Fi(x) by experiments, one can extract PDF at a certain Q 2 value. There can be two ways 
to compute PDF based on nonperturbative QCD techniques. By defining the interaction vertices in the handbag 
diagram, i.e. the qqir and qqj vertices, from effective models manifesting the spontaneous breakdown of chiral 
symmetry (SBCS), such as the NJL model [12], light-front formalism [33] and so on, one can compute directly the 
handbag diagram in the left of Figure [IJ using the optical theorem and the Bjorken limit Q 2 — »> oo, then determine 
F\{x) and PDF via Eqs. and ^ [12] [33]. It turned out that this method for computing PDF is equivalent to the 
operator-product-expansion (OPE) method, which will be explained below [34] . 

The other way is to extract the soft part for PDF from the hadron tensor W^ v in terms of the OPE technique 
considering the factorization theorem [15j [16]. Due to the OPE, the handbag diagram can be separated by the 
convolution of the hard and soft parts, and the separated diagrams are depicted in the right of Figure [T] The 
hard one is represented by Wilson coefficients, whereas the soft one by the matrix elements of local operators in a 
nonperturbative manner. In the present work, we employ this method with help of the gauge-invariant NLChQM. 
Now, we want to focus on the effective local operators for OPE, categorized by its twist and denoted by the gray blob 
of the soft diagram in Figure [l] Since higher twists are suppressed at high Q 2 , we will take into account only the 
twist-2 operator for PDF. The relevant matrix element of the local operator sandwiched by the PS-meson states can 
be related to the m-th moments of PDF as follows [15] : 



•ra+l 



-(<i>(p)\QfW(n ■ Drq f (0)\4>(p)) = / dxx m U(x) = (x^ 



/' 
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(3) 



where and n M denote the covariant derivative and light-like vector, which satisfies n 2 = 0, respectively. Note that 
the soft part is characterized by a renormalization scale /i, at which nonperturbative natures are determined. Hence, 
we will take \i implicitly in PDF hereafter. The 0-th moment with m = becomes for instance: 



-<#P)l9/(0)fa/(0)|<Kp)> 



/' 

^0 



dxf^x) = 1. 



(4) 



Note that Eq. (4| must be unity to satisfy the normalization condition for VQDF (or PDF). Readers can refer to 
Eqs. ([6| and (7) for the relations between various PDFs and VQDFs for tt and K. Alternatively, PDF can be also 
defined by pjlklES]: 



(5) 
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One can easily show that Eq. ([5| is equivalent to Eq. Q. The displacement was assigned by r] in Eq. Q. It can 
be understood by comparing the handbag and soft diagrams in Figure [I] the loop momentum k must satisfy the 
following condition, xp-n = k-n, after the factorization. This condition can be expressed in terms of a delta function 
in the loop integral over /c, i.e. S(k • n — xp • n), which corresponds to the effective composite local operator. At the 
same time, the momentum fraction x is also defined by x = (k • n)/{p • n). Here, the light-like vector n picks up the 
momentum (spatial) component of a vector, i.e. n • v = v + in the light-cone coordinate. In our theoretical framework, 
the following relations are satisfied for the positively charged pion and kaon: 

u<k(x) - d n (x), u^(x) = d n (x) — s 7T (x) = Snix) — 0> uk(x) = dx(x) = dx(x) = sk(x) = 0. (6) 

Here, the first relation in Eq. |6| comes from the isospin symmetry, m u = md, assumed in the present work. According 
to Eq. (J6|, we have the following relations between PDFs and VQDFs: 

v%(x) = Un(x), d%(x) = -d n (x) = -u n (x), sXix) = 0, v%(x) = u K {x), d%(x) = 0, s^(x) = -s K (x), (7) 

Now, we are in a position to compute the soft diagram representing PDF using an effective model. For that end, 
we employ NLChQM [25j [26] , being motivated by the instanton physics, which is properly defined in Euclidean 
space [33 [38]. The dilute instanton-liquid model is characterized by the nontrivial interactions between the quarks 
and (anti)instantons, which represents the nonperturbative QCD configuration, being characterized by the average 
inter- (ant i) instanton distance R ~ 1 fm and (anti) instanton size p ~ 1/3 fm. According to the nontrivial interactions, 
the quarks acquire their effective masses, resulting in the spontaneous breakdown of chiral symmetry (SBCS). The 
effective action of the instanton model can be bosonized, resulting in an effective chiral model with the quark and 
PS-meson degrees of freedom [38 . This effective chiral model possesses very interesting features: As mentioned, 
the dynamically-generated effective quark mass becomes (300 ~ 400) MeV to satisfy phenomenology and becomes a 
decreasing function of the momentum transfer. Hence, it plays the role of a UV regulator by construction, and the 
interaction strengths between the quarks and PS mesons becomes nonlocal, i.e. momentum dependent. Since we are 
interested in the physical quantities in Minkowski space, we perform the Wick rotation for the instanton model, i.e. 
t — >• it, then, obtain NLChQM, which resembles the nonlocal NJL model in many aspects. Although there is no firm 
theoretical proof for the validity of this analytic continuation between the instanton model and NLChQM, from a 
practical point of view, we employ NLChQM in the present work keeping the typical features of the instanton model. 

Taking into account all the ingredients discussed above, the effective chiral action (EChA) of NLChQM can be 
written in a neat form: 



S e ff [0,ra/;/i] = -iSpln 



■m f -y/M f (ip)U 5y /M f (ip) 



(8) 



where Sp stands for the functional trace, Tr c j n (x\ • • • in which (c, /, 7) denote the (color, flavor, Lorentz) indices, 
respectively, whereas rhf indicates the current-quark mass matrix, diag(ra n , m^, m s ). Throughout the present work, 
we will make use of the following numerical values for each mass: (m U: md,m s ) = (5,5,100) MeV [39], taking into 
account isospin symmetry and explicit SU(3) flavor- symmetry breaking. Mf is assigned for the nonlocal (momentum- 
dependent) effective quark mass. The nonlinear expression for the PS-meson fields, U§ reads: 



TT \ ^ X ' <M \ A 
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(9) 



Here, and A are the weak-decay constant for the PS meson as a normalization constant and the Gell-Mann 
matrix. One can easily see from Eq. ([8| that the nonlinear PS-meson fields and the quarks interact with the nonlocal 
(derivative) strength oc Mf. For instance, the effective Lagrangian density for the qqcj) vertex can be obtained from 
EChA as follows: 



-•nonlocal 



-q 



M f {d)>ft(\><l>)y/M f {d) 



q- (10) 



If we turn off the momentum dependence in Mj, it becomes a positive constant, and one can easily obtain the usual 
local PS-type effective Lagrangian density, 
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where stands for a massive constant. Since EChA in Eq. ([8| contains derivatives for the quark kinetic part as well 
as the quark effective mass, as suggested in Refs. [25], the gauge- invariant EChA can be easily obtained by imposing 



the minimal substitution in Eq. 



da 



iV^, in which V indicates a local vector field: 



Seff^TO/, V M ;/z] = -iSpln 



ip-rhf- JM f (ip)U 5 JM f (ip) 



(12) 



In this way, one ensures the gauge invar iance of relevant physical quantities, extracted from Eq. ( |T2| ). We also note 
that this gauging procedure is similar to that given in Ref. [15 , in which the single-instanton model was employed. 
Here is one theoretical caveat: When our writing the gauge- invariant EChA as in Eq. (12), we assumed the gauge 
connection in a straight-line path, which is most convenient and practical [40] . 

Considering the effective local vertex for m = in Eq. (|3| , the matrix element in the right-hand-side of Eq. ([5| for 
PDF can be evaluated by the three-point functional derivative with respect to the (j) and V in Eq. ( 12 ) with the delta 
function, 



S^(x)5^(y)5V^0) 



(13) 



where the superscripts (a, /?) denote the isospin indices for the PS mesons. Simultaneously, we expand the nonlinear 
PS-meson field U$ in EChA up to 0(4> 2 ), since we are interested in the two PS meson fields for the initial and final 
states as in Eq. After performing these procedures, the analytical result for PDF via NLChQM reads: 




cTfc 



5{k a -n — xp- n)Tr 7 



M b "f 5 ^/MaSaflSa ^ M al5 ^ M b S b 

M 675 ( v/Ma -n)S a v/Ma75 \ZM~bS b 



(14) 



The subscripts (a, b) indicate two different flavors inside the PS meson with the momenta (k a ,kb) = (/c, k — p) as 
defined in Figure [I] We note that the delta function, S(k a • n — xp • n) has been convoluted in the integral over k as 
mentioned above. It is also emphasized that the second and third terms in the square bracket in the right-hand-side of 
Eq. (14) exist only when we consider the momentum-dependent effective quark mass, since those terms are obtained 



from the functional derivative of the gauge-invariant EChA in Eq. ( 12 ) with respect to which is a special feature 



of the present nonlocal interaction model. In other words, these derivative or nonlocal contributions ensure the gauge 
invariance for PDF by construction [25]. The quark propagator for the flavor a is denoted by 



S a 



[h + (m„ + M f )} 



Ih + Mf] 



kl - (m a + M a y + ie kl-Mi + ie' 



(15) 



In the above equation, we introduced a notation M a and indicated the Feynman e explicitly. The relevant mass 
functions in Eq. ( 14 ) read 



M a = M 



k 2 a -/i 2 



a (kl - /i 2 + ie) ' 



(16) 



Note that we employ the Lorentzian-type structure function for the effective quark mass, as motivated by the instanton 
physics [38] as well as employed in Refs. [25-27 . Here, Mo indicates the constituent-quark mass at zero virtuality 
which is the model parameter to be determined phenomenologically in the next Section. Evaluating the trace over 
each spin index and employing the light-cone coordinate vector manipulations 



k ■ n — k + = xp + , k 2 = k + k 



fc T , p 



k • p = 



p + k + k + p 



71+ 



(17) 



we are left with the contour and polar integrals for k and fc^, after the delta function integral over k + in Eq. (14): 

iN c f dk~d 2 kT r ^ /,- 7 2\ t~ /,_ ,2n ~r~ n- 7 2\i / -\ 
"4F 2 J (2tt) 3 I ' ,fcr)+^NL, a (fc ,fcr)+^NL,b(* , k T )\ + (x o x), 



U( x ) 



(18) 



where the relevant functions *7-L,NL,(a,b) m the square bracket are explicitly defined in Appendix. We have also used 
a notation x = (1 — x) here. The integral in Eq. ( [18] ), however, is not simple in comparison to those of usual local 
models. Since M a ^ in Eq. (16), which can have poles in the denominator during the integrals over appear in the 
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denominator as well as the numerator in the integrand of Eq. (18) simultaneously. For instance, we have to solve 
the septic equation of k~ in the present work to find the poles in performing a contour integral over k~ with an 
appropriate cut. The choice of the cut is discussed in detail in Ref. [26], and we follow their cut scheme. It is worth 
mentioning that the numbers of the poles appearing in the calculation of PDF relates to the power of Mf in Eq. (16). 
For instance, if we take Mf = M [- • • ] lor3 , it is necessary to solve a quintic or nonic equation to find the poles in the 
present nonlocal model. As shown in Ref. [25], the change of the power can make effects on the shape of the curves. 
We will, however, only consider the case with Eq. (16) in the present work for simplicity, since the changes of the 
shapes due to the different power appear only around the end points, x = and x = 1, [25] and the changes can be 
absorbed in the DGLAP evolution qualitatively. After the contour integral over k~ , the polar integrals over kr can 
be done easily numerically, and one is led to the final result for PDF (or VQDF). 



III. NUMERICAL RESULTS AND DISCUSSIONS 



In this Section, we provide the numerical results with relevant discussions. First, we would like to explain how to 
determine the model parameters, saying, the model renormalization scale \i and the constituent-quark mass at zero 
virtuality Mq in Eq. (16). For this purpose, we make use of the PDF normalization condition in Eq. Q and the 



empirical values for the PS-meson weak-decay constants, F^,k = (93.2, 113.4) MeV. Note that F^ appears in Eq. (fl8| 
as a normalization constant in the denominator. Considering the phenomenological scale for the hadrons aroundl 
GeV, first, we try \i = (0.8 ~ 1.2) GeV to satisfy the normalization condition with F n = 92.3 MeV as an input for the 
pion PDF, by varying the M value. Corresponding numerical results for u^{x) are shown in the left panel of Fig. [2] 
As for the various renormalization scales, Mo varies from 260 MeV to 360 MeV to satisfy the normalization condition. 
Moreover, depending on the parameter set of (/i,Mo), the shape of PDF also changes significantly. We observe a 
tendency that as \i increases (and vice versa for Mo), PDF gets broadened. In order to choose one parameter set out of 
the five as given in the left panel of Fig.JiJ we also consider the kaon case similarly. We compute Fk = i^? mputed 5 using 
those parameter sets and Eq. Q, then compare them with its empirical value. One can choose the best parameter 
set, which minimize the quantity, defined as 

A zt 1 I 7-iComputed TTiempiricai i / in \ 

= \r K ~ r K V (19) 

We list the numerical results for Fk in Table [ij As shown in the table, the parameter set (/i,M ) = (1, 0.3) GeV 
presents the smallest deviation ~ 7% from the empirical value of Fk- Although we note that the deviation can be 
reduced by tuning the current-quark masses (m n ^,m s ), being different from (5, 100) MeV, we do not perform those 
tunings, since they only produce qualitatively negligible changes in the numerical results. Thus, we will choose this 
parameter set for all the numerical calculations hereafter. Using the parameter set chosen, we show the numerical 
results for u n (x) for the total (solid), local (dot), and nonlocal (dash) contributions separately in the right panel of 
Fig. [2] We observe that the local contribution has a peak x = 0.5, whereas the nonlocal one show bumps around 
x = 0.3 and 0.7. By integrating the local and nonlocal contributions over x, we have about 0.76 and 0.24, respectively. 
The ratio 0.24/0.76 ~ 1/3 is a typical value for each contributions in the NLChQM calculations, when the gauge 
invar iance is taken into account explicitly. 

In the left panel of Fig.|3j we draw the numerical results for u^{x) (solid), uk{x) (dot), and sk{%) (dash). Due to the 
far heavier mass of the strange-quark mass ~ 100 MeV, the curves for the kaon present obviously asymmetric shapes, 
signaling the explicit flavor-SU(3)-symmetry breaking, while that for the pion is clearly symmetric with respect to x 
according to the isospin symmetry. Since, as in many literatures to investigate PDFs, it is convenient to use the form 
xfcfyfx) rather than /^(x) [T6J |22j |23l [41] , and the numerical results for those functions are given in the right panel of 
Fig. [3] We note that all the results are qualitatively compatible with other theoretical results [T2 | [T5HT7] . Moreover, 
the shapes of the curves are different from the previous works [22j [23] , indicating that the explicit considerations on 
the momentum-dependent quark masses in the calculations are crucial, in addition to the gauge invariance. 



fi [GeV] 
M [MeV] 


0.8 
360 


0.9 
314 


1.0 1.1 
300 275 


1.2 

260 


^computed [MeV] 

AF K [MeV] 
Deviation [%] 


126.1 
12.7 
11.2 


122.3 
9.0 
7.9 


121.7 122.2 
8.3 8.8 
7.3 7.8 


123.0 
9.6 
8.5 



TABLE I: Computed Fk with various par ameter sets (/i,Mo), which are determined by Eq. Q and the input = 93.2 MeV. 
We also list the values for AFk in Eq. |l9| ) and deviations from its empirical value. We will use the third parameter set, written 
in a boldface, for all the numerical results. 
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FIG. 2: (Color online) Left: PDF for the pion, u^{x) for various renormalization scale fi and constituent-quark mass at zero 
virtuality Mo, satisfying the PDF normalization condition Eq. Q with — 93.2 MeV. Right: u-n(x) is shown separately for 
the total (solid), local (dot), and nonlocal (dash) contributions for (/z, Mo) = (1,0.3) GeV. 

i 1 1.4 i 1 




FIG. 3: (Color online) Left: Valance-quark parton distribution functions (VQDF), u 7r (x) (soild), uk(x) (dot), and sk(x) (dash) 
from Eq. (18). The vertical lines indicate the position for the momentum fraction x = 0.5. Here, we use (/i, Mo, F^, Fk) ~ 
(1000,300,93.2, 121.7) MeV and in order to satisfy the normalization condition for PDFs in Eq. Q. Right: Those multiplied 
by x, i.e. xu^^x), xuk(x), and xsk(x) are also given in the same manner. 



The parameterization of PDF in the following analytic form is very useful for various applications in data analyses 
and high-Q 2 evolutions: 

xU(x) = a^x b *(l-x) c ^ (20) 

where the fitting parameters (a, 6, c)^ are certain positive-real values. Using the numerical results shown in the right 
panel of Fig. |3j those parameters can be fitted and are listed in Table. [XX] for the low-renormalization scale fi = 1 
GeV. Using these parameterized PDFs in Eq. (20) and Table [n] as inputs, we perform the high-Q 2 evolution, using 



the DGLAP equation, to compare our results with the empirical data and theoretical estimations in Refs. [29Hc 
The fortran code QCDNUM is used to this end [42j 03]. In the left panel of Fig. ^ we show xu£(x), evolved to 
Q 2 = 27 GeV 2 using the LO (thin) and NLO (thick) DGLAP evolutions. Here, the empirical data for the pion 
are taken from the muon-pair production experiment by 252 GeV pions on tungsten target [29 . Since it is rather 
uncertain to choose the initial momentum scale Qo in the DGLAP evolution, we try three different values for it as 
Ql = (0.15, 0.20, 0.25) GeV 2 . Note that these values are related to the momentum scale about (400 ~ 500) MeV, 
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FIG. 4: (Color online) Left: Valance w-quark distribution function multiplied by x, i«„ 0) evolved to Q 2 = 27GeV 2 for various 
initial values Ql = (0.15, 0.20, 0.25) GeV 2 in (solid, dot, dash) lines, using the LO (thin) and NLO (thick) DGLAP evolutions. 
Empirical data are taken from the muon-pair production experiment by 252 GeV pions on tungsten target [29]. Right: The 
same for xu K {x) with the NLO DGLAP evolution. 



which are also compatible to typical nonperturbative scales A = (0.5 ~ 0.6) GeV in the instanton model [37J |38] and 
A w 0.4 GeV in the generic NJL model [44]. Again, these values are very close to those used in Refs. [12] HZ] [45] 
for the same purpose. As shown in the left panel of Fig. [4j the numerical results reproduce the empirical data 
qualitatively well for the NLO DGLAP evolution, whereas the LO DGLAP evolution results in general show larger 
than the empirical data. Since the NLO DGALP results are much more compatible to the data, we will present the 
numerical results only for the NLO ones hereafter. Similarly, we present the numerical results for xu^(x), evolved 
to Q 2 — 27 GeV 2 , in the right panel of Fig. [ZJ although there have been no empirical data to be compared to. In 
comparison to xv%(x), the peak positions are shifted to the lower x region, due to the explicit flavor- SU (3) -symmetry 
breaking, as understood in the right panel of Fig. [3] Although there are strength differences in the VQDF curves for 
each Qq values, the curve shapes are qualitatively similar to each other. 

Now, we present the numerical results for the ratio of two different PDFs, i.e. u n (x)/u K (x) at Q 2 = 27 GeV 2 for 
the three different initial values Q 2 , = (1.5,2.0.2.5) GeV 2 (solid, dot, dash) in the left panel of Fig. [5] They can be 
compared with the data from the 150 GeV incident-beam experiment for (K~ , tt~) + nucleus — ^ /i + \i~ X [46]. We also 
depict the fitted curve for the data: uk(x) /u n (x) = 1.1(1 — x) 0,22 (long dash) [47]. Although the numerical results 
match with the data qualitatively well in the region x > 0.5, we observe overshoots for the smaller x by (20 ~ 30)%. 
Moreover, there are visible deviations between the numerical results and the fitted curve of Ref. [47 . We note that 
the curve shapes of our results resemble that from the reduced Bethe-Salpeter equation (BSE) vertex calculation [T7] . 
although the strengths are slightly different. Except for the region x < 0.2, the result given in Ref. [12 is also similar 
to ours in shape. In the right panel of Fig.[5j the numerical results for sk(x)/uk(x) are given in the same manner with 
the left one. Interestingly, the curve behavior is quite different from that for uk{x)/u^{x). Again, this discrepancy 
can be understood by the explicit flavor-SU(3)-symmetry breaking effects. As observed by comparing the curves in 
Figs. and ffl fY(x) in the vicinity of x = 1 is almost unaffected by the DGLAP evolution. Hence, as discussed in 
Refs. [12] [T7]T f¥(x — » 1) can be considered to represent the relatively pure nonperturbative QCD contributions. In 





CLcf) 


bcf> 


C0 




7.60 


2.16 


1.13 


UK 


18.90 


2.43 


2.00 


SK 


35.67 


3.47 


1.83 



TABLE II: Parameters (a, b, c)^ for the pion and kaon PDFs multiplied by x in Eq. (20). 
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FIG. 5: (Color online) Right: Ratio uk/u-k at Q 2 = 27GeV 2 as a function of x for various initial values Qq — 
(0.15, 0.20, 0.25) GeV 2 in (solid, dot, dash) lines. We also show the fitted curve in Ref. 47 (long dash). Empirical data 
taken from the 150 GeV incident-beam experiment for (K~ ,tt~) + nucleus — > /i + ji~ X |3B]. Right: The same for sk{x)/uk{x)\. 



the previous work [23] , the ratio of uk/u^ was approximated for the leading local contribution as 



u K {x) 



u n (x) 



F 2 k 



1 



Mo 



0.4. 



where we have used relevant quantities employed and computed in the present work, i.e. 

(F„, F K , M , m u , m s ) = (93.2, 121.7, 300, 5, 100) MeV. 



(21) 



(22) 



Considering that the nonlocal contribution produces (20 ~ 30)% additional strength, we have approximately uk(x — >• 
1)/utt(x 1) = (0.52 ~ 0.53), which is comparable with the actual value 0.45. It is worth mentioning that different 
forms for the ratio were suggested as (M u /M s ) 2 from the NJL calculation [12] and (F n / 'F K )(M U /M S ) 4 from the BSE 
calculation [T7] . 

Finally, we would like to discuss the moments of the pion and kaon VQDFs as defined in Eq. Q. For brevity, we 
will take into account those only for the u quark. They are very useful to analyze the data and determine unknown 
parameters in models and the DGLAP evolutions. In Refs. 
given as follows: 



, the first two moments for the pion, (x 72-1 ' 2 )^ are 



2(x) 7T 
2(x) 7T 
2(x) 7r 



0.55, 2(x 2 ) 7T « 0.18, 



0.40 ±0.02, 
0.46 ±0.07, 



2<* 2 > 
2(x 2 ) 



(fitl) : Q 2 =4 GeV 2 [31], 
= 0.16 ±0.01, : Q 2 =4GeV 2 [30], 
= 0.18 ±0.05, : Q 2 =49GeV 2 [32]. 



(23) 



The values in Ref. [31] were obtained by the NLO analysis (fitl) of the data for the Drell-Yan process tt~ N — >• /i + \i~ X 
of the experiments by E-615 [29] and NA10 [48] collaborations. Ref. [30] also analyzed the data of NA10 collaboration. 
On the contrary, the moments in Ref. [32] were obtained by the lattice QCD simulation with the Wilson fermions 
in the quenched approximation for the pion structure functions. From our numerical results, we have the followings, 
using the curves depicted in Fig. [4j 



2(x) 7T = (0.37 - 0.46), 2(x 2 



(0.14-0.19) : Q 2 =27 GeV 2 



(24) 



for Qq = (0.15,0.20,0.25) GeV 2 as shown in the left panel of Fig. [HJin the (square, circle, triangle), respectively. We 
note that, although the Q 2 values are different from other estimations, these values in Eq. (24) are well compatible 
to those given in Eq. ( [23] ) . The schematic comparison of these moments from different works are given in the left 
panel of Fig. [6] For clearance, we present the data of Refs. [30] (diamond) and [32] (nabla), shifted by —0.5 and ±0.5 
for n, respectively. Since, in Ref. [31], the authors provided the parameterized VQDF as u%(x) = 0.077 x x _0,85 (l — 
x) 1,75 (l ± 89 Ax 2 ) at Q 2 = 4 GeV 2 , we can compute other higher moments beyond the first two as in the left panel of 
Fig. [6] (rhombus) . As shown there, the present numerical results for the moments are in qualitatively good agreement 
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FIG. 6: (Color online) Left: n-th moments for the pion VQDFs for Ql = (0.15 - 0.25) GeV 2 at Q 2 = 27 GeV 2 . The data 
are taken from Aicher et al. [31 , Martinelli et al. [32], and Sutton et al. [30] . For clearance, we shift the data of Ref. [32] and 
Ref. [30] by -0.5 and +0.5 for n, respectively. Right: The same for the kaon. For the numerics, see Table III 



with other analyses and theories. In the right panel of Fig. |6j we also present the moments for the kaon VQDF, 
as functions of n. Overall behavior of the curve is very similar to that for the pion, whereas visible differences are 
observed for small n values. This tendency can be understood clearly by seeing Table |III[ in which we list all the 
moments for the pion and kaon for those Qo values at Q 2 = 27 GeV 2 . From the numerics given in the table, we see 
that the higher moments for n > 2 are very small and almost identical for the two mesons. As a consequence, we 
have the following observation: 

/ n=l,2\ /_n>2\ 

\ x )J£ ~ I I UL ~ 1 (25) 



IV. SUMMARY AND CONCLUSION 



We have investigated the pion and kaon PDFs employing the OPE technique for the handbag diagram with the 
twist-2 operator and the gauge- invariant NLChQM, which conserves the vector currents. All the model parameters 
were determined by the empirical inputs and model-independent constraint, i.e. PS-meson weak-decay constants and 
the normalization condition for VQDF. Once all the parameters fixed, we performed the high-Q 2 evolution by the LO 
and/or NLO DGLAP equation for the different initial momentum scales. In what follows, we summarize important 
observations in the present work: 

• In order to satisfy the gauge invariance, we include the nonlocal contribution to PDF in addition to the leading 
local one. The nonlocal contribution provides about 30% strength for PDF, as usual in the similar approaches. 
Appropriate consideration for the momentum-dependence of the effective quark mass changes the shapes of 
PDF, in comparison to the simplified calculations as in the previous work. 





Ql [GeV 2 ] 


n 1 


23456789 10 




0.15 
0.20 
0.25 


0.184 
0.214 
0.230 


0.068 0.033 0.018 0.011 0.007 0.005 0.004 0.003 0.002 
0.087 0.044 0.026 0.016 0.011 0.008 0.006 0.004 0.003 
0.097 0.050 0.030 0.019 0.013 0.009 0.007 0.005 0.004 




0.15 
0.20 
0.25 


0.207 
0.242 
0.261 


0.072 0.032 0.017 0.010 0.006 0.004 0.003 0.002 0.002 
0.091 0.043 0.024 0.014 0.009 0.006 0.004 0.003 0.002 
0.102 0.050 0.028 0.017 0.011 0.007 0.005 0.004 0.003 



TABLE III: n-th moments for the pion and kaon VQDF for Ql = (0.15 - 0.25) GeV 2 at Q 2 = 27 GeV 2 . 
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• Taking = 93.2 MeV as an input, we determine the renormalization scale \i — 1 GeV as well as the constituent- 
quark mass at zero virtuality Mq = 300 MeV, together with the information Fk = 113.4 MeV. These values 
give Fk — 121.7 MeV, which is slightly larger than its empirical value by about 7%. Computed PDFs for the 
pion and kaon at the low-renormalization scale are compatible qualitatively with other theoretical results, such 
as the single-instanton model, BSE approach, NJL model, and so on. 

• After performing the DGALP evolution up to Q 2 = 27 GeV 2 for three different initial momentum scales Qq = 
(0.15, 0.20, 0.25) GeV 2 , we compare VQDFs for the pion and kaon with the empirical data as well as theoretical 
calculations, resulting in qualitatively good agreement with them. However, we still observe sizable differences 
with the empirical data for some specific x regions, especially as shown in uk(x)/u 7T (x). 

• It turns out that the moments for the pion and kaon VQDFs are well compatible with experimental and 
theoretical estimations. We observe that the differences between the higher moments for n > 2 are almost 
negligible, whereas the first two moments, i.e. n = (1,2), show a tendency that (x n=1 ' 2 )x/(x n=1 ' 2 ) 7r « 1.1. 

Since we are now equipped with PDF in the gauge- invariant manner, it must be interesting to apply this PDF to 
the investigations for other physical quantities. For instance, we can use the present result to compute the nucleon 
PDF together with the PS-meson FF, together with an ansatz for the bare nucleon PDF as studied in Refs. [49U5T] . 
Related works are under progress and appear elsewhere. 
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Appendix 

The relevant functions in Eq. (fl8|) are defined as follows: 



4p+ri 2 D$ D A a D s h [xk 2 T + x 2 k~p+] + 2xD\ {p\m h + rj) (D A a m a + rj) + xD\ (D*m b + rj) 



Ap+rf 



D*(£ a - ak~) + 2m aV Di + ri 2 ] a [D*[£ h - f3k~ + 6)] + 2m b r 1 Dl + rj 2 Y h 
rj (D%m b + rj) + D^[m a r] + B\(2k\ + xk~p+ + m a m b + xm 2 ^) 



[D*(£ a - ak~) + 2m a rjDi + rj\ [D*[£ b - j3k~ + 6)] + 2m b rjDt + rj 2 ] b 



4p+rj 2 



rj (D b m b + rj) + D^[m a r] + D\(2k\ + xk p + + m a m b + xmV)] 



[D*(£ a - ak~) + 2m a rjDi + V \ [D*[£ b - j3k~ + 6)] + 2m 6 ^ + 



7 — afc~ 



7 - f)k- + 5 



where represents the PS-meson mass and the various notations read: 



— x, 2 — x, 1 — 2x), 

(a, (3, 7, 8, t], £ a , £ b ) = (xp + , -xp + , k^ + n 2 , -xm\, Mq/i 4 , k^+m 2 a , k% + m\), 
(D 2 a , D 2 b ) = (7-afc", 1 -pk~+5). 



(26) 



(27) 



In deriving above equations, we have assumed that p 2 = m 2 = p + p~ — p\ ~ P + P~ , considering that the transverse 
momentum for the meson is much smaller than the longitudinal ones. i.e. pt ~ . 

In order to perform a contour integral for Eq. (|26| ) over k~ , we pick the poles from the denominator. We have 
['■']« = Si=o a i(k~Y m the denominator of Eq. (I26F7 The coefficients a$ are written as 



a 5 = a 5 , a 4 = -a 4 (4 7 + £ a ), a 3 = a 3 (6 7 2 + 4 7 £ a ), a 2 = -a 2 (4 7 3 + 6 7 2 £ a + 2ra a 7?), 
ai = a 1 (7 4 + 4 7 3 £ a +4ra a 777), a = -a°(7 4 ^a + 2m a j 2 r] + r? 2 ). 



(28) 



The poles can be obtained by solving this quintic equation with respect to k numerically. 
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